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ABSTRACT 


This  report  provides  estimates  of  the  positional  accuracy  for 
celestial  objects  which  we  have  obtained  at  the  Experimental  Test 
Site  of  the  Ground-based  Electro-optical  Deep  Space  Surveillance 
(GEODSS)  program.  Optical  observations  made  at  the  ETS  have  been 
compared  with  radar  measurements  performed  by  the  Millstone  Hill 
Radar.  A brief  description  of  the  equipment  and  theoretical 
basis  of  the  optical  data  reduction  is  included  too.  It  appears 
that  optical  observations  of  artificial  satellites,  with  real 
time  reduction,  are  accurate  to  =5". 
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I . INTRODUCTION 

For  twenty  years  the  Baker-Nunn  camera  system  has  admirably  served  as 
our  principal  source  of  optical  obsei-vations  of  artificial  satellites. 

The  United  States  Air  Force  and  Lincoln  Laboratory  have  recently  completed 
the  design,  construction,  and  testing  of  the  Experimental  Test  Site  ( /ITS) 
of  the  Ground-based  Electro -Optical  Deep  Space  Surveillance  (GEODSS)  program. 
GEODSS  is  now  currently  being  deployed  world-wide  by  the  TRW  Company.  When 
fully  operational,  it  will  supplant  the  Baker-Nunn  photographic  camera  system. 
Moreover,  because  of  various  technological  improvements,  its  limiting 
magnitude  Gd;16m),  search  rate,  and  observing  rate  will  far  exceed  those 
of  the  Baker-Nunn  camera  network.  Hence,  we  can  hope  that  some  telescope 
time  can  be  devoted  to  purely  astronomical  research  (e.g.,  minor  planet 
observations,  supemovae  searches,  deep,  multicolor,  and  very  rapid 
catalogue  construction  to«l”,  variable  star  searches,  meteor  counts,  etc.). 
For  these  reasons  we  believe  that  the  positional  accuracy  of  GEODSS' s 
observations  will  be  of  interest  to  a large  audience.  Hence,  this  publica- 
tion. 

Below  we  briefly  describe  the  hardware  and  computer  facilities  at 
the  ETS.  We  then  touch  on  some  purely  astrometrical  matters  and  our 
handling  of  the  SAOC  (which  serves  as  our  reference  star  catalogue). 

Next  we  discuss  the  theoretical  basis  for  the  Precision  Local  Calibration 
(PLC)  procedure  currently  employed  by  tiie  ETS.  Ibis  procedure  rests 
directly  on  the  reduction  techniques  of  classical  photographic  astrometry. 
After  presenting  PLC  we  consider  its  expected  accuracy  and  some  simulation 
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testing  we  have  nerformed.  Observations  made  in  May  of  1978  on  the 
brighter  PK4  stars  and  eight  artificial  satellites  have  been  used  to 
experimentally  verify  the  analytical  and  numerical  work.  Lastly,  we 
describe  a simpler,  less  accurate  Single  Star  Calibration  procedure 
which  has  been  tested  in  the  same  way.  Ultimately,  we  expect  to  use  the 
plate  overlap  technique  in  real  time  (Taff^). 

Briefly  the  results  are  as  follows:  The  positional  accuracy  for 
the  stars  is  173  +_  170  and  -5"  + 3"  for  the  satellites.  Secondly, 
relative  to  the  FK4  itself  for  the  stars  or  the  Millstone  Hill  Radar  for 
the  artificial  satellites,  there  is  no  systematic  bias  in  right  ascension 
or  declination.  Thirdly,  there  is  no  measurable  color  term,  magnitude 
term,  or  systematic  star-satellite  term.  In  addition  to  these  results, 
we  wish  to  stress  that  the  entire  procedure  (e.g.,  PTC  or  SSC)  occurs  in 
real  time,  under  computer  control,  with  an  electronically  generated 
finding  chart,  and  requires  approximately  ll’s  (21s  for  SSC)  to  complete. 
Furthermore,  our  theoretical  computations  and  telescope  alignment  clearly 
indicate  that  we  are  only  limited  by  the  hardware  from  routinely  splitting 
the  seeing  disc. 


II.  EQUIPMENT 


'Hie  telescope  is  a 31"  Boiler  and  Chivens  f/5  of  Ritchey-Chretien  design. 

At  the  Cassegrain  focus  the  field  of  view  is  1?16.  A zoom  feature  of  the 

camera  reduces  this  to  34 ! 8 (2:1).  The  low  light  level,  high  gain  camera 
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was  built  by  Westinghouse  following  a design  of  R.  Weber. ~ It  consists  of 
a Varian  VLI-116  image  intensifier  coupled  via  fiber  optics  to  a Westing- 
house  WX-32719  camera  tube.  The  first  photosurface  is  an  80mm  diameter 
S-20  type.  The  final  target  is  a 32mm  diameter  silicon  diode  array. 

Standard  U.S.A.  television  frame  rates  and  horizontal  scan  line  densities 

were  used  for  the  work  discussed  herein.  The  camera  target  has  450  resels/line 

and  the  right  ascension:  declination:  field  of  view  ratio  is  4:3:5.  In 

the  zoom  configuration,  the  resolution  element  size  is  3'.'7.  The  50" 

response  area  has  a diameter  of  56p  on  the  camera  target.  The  high 

(5  x 10^)  pre-readout  gain  ensures  background  noise  limited  performance 

with  uncooled  electronics. 

The  control  of  the  telescope,  the  storage  of  the  data  files,  and  the 
data  reduction  are  all  handled  by  a MODCOMP  IV  computer.  All  the  data  files 
(star  catalogues,  Besselian  day  numbers,  etc.)  are  stored  using  a 32  bit 
word  for  right  ascension  and  declination  (O'.'IO  in  each  coordinate).  All 
arithmetic  is  performed  using  a 48  bit  word  length. 

The  operator  views  the  telescope  display  on  a television  monitor  with 
a set  of  rectangular  gridlines  super impo s ed . While  all  gross  movements  of 
the  telescope  are  under  computer  control,  the  final  alignment  of  the  star 
(or  artificial  satellite)  relative  to  a fiducial  mark  in  the  grid  system 
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is  performed  using  a joystick  type  device.  This  control  allows  one  to 
select  one  of  2 = 256  different  horizontal  and  vertical  settings  and 

provides  a 675  positioning  uncertainty  independent  of  the  camera's  resol u 
tion  element. 

A statistically  independent  combination  of  tiie  telescope's  shaft 
encoder's  resolution  (0762),  the  atmospheric  seeing  disk  (=1"),  the  camera's 
resolution  cell,  and  the  joystick's  resolution  cell  imply  that  the  standard 
deviation  of  a single  measurement  is  - 5”.  Experiments  designed  to  test 
this  result  have  indicated  that  this  figure  is  probably  closer  to  4". 

The  observatory  is  located  on  the  White  Sands  Missile  Range  near 
Socorro,  New  Mexico.  Additional  information  about  the  observatory  is 
contained  in  reference  2. 
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III.  POSITION  REDUCTIONS 


The  source  for  positions  and  proper  motions  was  the  Smithsonian 

Astrophysical  Observatory  Star  Catalog  (Staff  of  the  SAO,  1966).  These 

1950.0  values  were  updated  for  precession  and  proper  motion  to  1978.0  using 
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the  power  series  method  (Newcomb,  Fricke  and  Kopff  , and  Taff  ).  Foreshorten 

ing  terms  were  not  included.  The  accuracy  is  at  least  OT’OOl  in  right 

ascension  and  O'.'Ol  in  declination  for  polar  distances  larger  than  9° 

(exclusive  of  the  foreshortening  terms). 

The  reduction  of  the  reference  stars  from  catalogue  mean  place  to 

topocentric  mean  place  was  performed  using  the  Besselian  day  numbers  A-E, 

J,  and  J'  as  well  as  including  diurnal  aberration.  The  astronomical 

6 

refraction  correction  followed  Smart  . The  annual  parallax  correction 
was  not  included.  A totally  differential  reduction  from  mean  place  to 
topocentric  place  has  been  implemented  in  which  the  reference  stars  arc 
reduced  relative  to  the  center  of  the  field.  This  allows  one  to  recover 
the  program  object's  mean  place  immediately  after  its  topocentric  place 
has  been  obtained  from  the  telescope/camera  model.  It  was  these  results 
that  were  directly  compared  (in  the  case  of  the  stars)  to  1978.0  values. 

The  reduction  of  the  artificial  satellites'  position  is  performed 
identically.  Both  parallactic  refraction  and  planetary  aberration  are 
small  (=2")  and  presently  ignored.  Complete  details  are  in  reference  5. 
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IV.  THi:  TELESCOPE/CAMERA  MODEL 


Each  field  contains  N reference  stars  and  the  program  object.  The 
topocentric  position  of  a reference  star  is  denoted  by  fa, 6)  and  the  topo- 
centric  position  of  the  program  object  is  denoted  by  (A, A).  The  corres- 
ponding values  as  measured  by  the  telescope  are  (ap  6y)  and  (ApA,-) . 

Ihe  model  relates  { (a , 6) } to  {(aT,5^)}.  then  a prediction  of  (A, A)  is 
made  from  (ApAp)  by  inverting  the  model.  It  was  assumed  that  a linear 
model  would  suffice  . The  equations  below  are  analagous  to  the  ones  used 
when  modeling  a photographic  plate  in  classical  photographic  astrometry. 
The  intermediary  device  of  standard  coordinates  was  not  used. 

Let  an  angular  bracket  denote  an  equally  weighted  average  over  the  N 
reference  stars  within  a field.  Furthermore,  for  each  of  the  reference 
stars  define 


ea  a ' “T»  e<5  6 " 5T  ’ 

and  e = a - <a>,  = 6 - <6>. 

a ’ 5 

The  corresponding  quantities  for  the  program  object  are, 

E = A - <a>  and  E.  = A - <6>. 
a o 

Ihe  full  linear  "plate  model"  is 


e = ac 

+ ber  + c, 

a a 

6 * 

e r = dc 
o a 

+ ee^  + f, 

where  a - f are  the  "plate  constants". 


(la) 

(lb) 

(2) 

(5a) 

(3b) 


* 

Both  the  telescope  and  the  television  camera  arc  being  replaced  so  a 
detailed  analysis  of  the  geometry  of  the  interface  has  been  postponed. 


b 


1 


The  plate  constants  are  determined  by  the  procedures  of  least  squares 
with  equal  weighting  for  each  component  of  the  model.  Once  this  is  done, 


the  topocentric  position  of  the  program  object  i~  recovered  by  solving 
Eqs.  (4); 

A - ,'Vp  = aE^  + bEfi  + c,  (4a) 

A - AT  = dE  + eEx  + f.  (4b) 

1 a o 

The  determinant  of  the  normal  equations  is  given  by  Ndet(I)  where 


and 


I = Ee  e , I ■ 
aa  a a ao 


Ee  er . 
a 6’ 


Eexer. 
o o 


Hence,  I is  the  "moment  of  inertia  tensor"  of  the  reference  star 
distribution.  Since  det(l)  vanishes  if  and  only  if  the  reference  stars 
lie  on  a straight  line,  as  long  as  the  problem  is  well  posed,  the  solution 
for  the  plate  constants  exists  and  is  unique.  Furthermore , if  denote 

the  standard  deviations  of  a single  measurement  in  right  ascension  and 
declination,  then  from  the  covariance  matrix  of  the  plate  constants, 

var(A)  = {aa2  (1-e)  2 + a26b2Hl  + 1/N  + E • I"1  • E)/D2,  (6a) 

var(A)  = (oa2d2+  c*.  2(l-ai  2}{1  + 1/N  + E 'I'1  • E}/D?,  (6b) 
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where 


— (Ea’V’ 

and  D is  the  determinant  of  Eqs.  (4), 

D = (1-a) (1-e)  - bd. 

From  Eqs.  (6)  we  can  conclude  that  the  optimum  position  for  the  program 
object  is  at  the  center  of  the  field  (remember  I is  a positive  definite 
quadratic  form).  Moreover,  when  the  program  object  is  so  placed,  the 
expected  accuracy  is  independent  of  the  areal  extent  of  the  reference 
stars  and  increasing  the  number  of  reference  stars  beyond  the  minimum 
value  of  4 offers  little  improvement  in  accuracy  (its  2.3 % for  N=4  -*  N=5 
and  11.81  for  N=4  N=°°) . While  further  rigorous  analytical  progress 

appears  to  be  difficult,  a heuristic  averaging  over  the  noise  leads  one 
to 


< var(A)  > - o 2 

[i  + i/N  + e • r1 

a. 

< var(A)  > - o^  2 

[1  + 1/N  + E • T"1 

However , 

E • I"1  • E = E2a/var(ea)  + F. 2-/var(c(f ) , (8a) 

so  that  if  E f 0 one  wants  the  maximum  areal  extent  for  the  best  accuracy. 
In  fact,  if  the  right  ascension  and  declination  coordinates  are  treated 
symmetrically,  then 

E • I'1  • E ~ kE  2/(NS)  (8b) 

where  S is  the  area  occupied  by  the  reference  stars  and  k is  related  to 
the  radius  of  gyration  of  the  distribution  of  reference  stars  about  the 
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Extensive  Monte  Carlo  simulations  using  a uniform  distribution  of 


reference  stars  and  treating  right  ascension  and  declination  symmetrically 
have  confirmed  Eqs.  (7,8).  Additional  simulations  were  conducted  using 
the  SAOC  itself.  The  field  centers  for  the  later  computations  were  at 
a = c/1(ol15)23*}5,  6 = -85°(5°)850.  While  the  SAOC  simulations  only  used 
a circular  field  of  0?75  radius,  the  Monte  Carlo  trials  used  rectangular 
fields  of  edge  size  0?25,  0?50,  1?00,  2?00,  and  4?00.  Foi  both  sets 
of  computations  c = a,;  = a with  a = 1725,  2750,  5700  or  10700  and  the 
values  of  N used  were  3(1)8.  See  reference  9 for  complete  details. 
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V. 


RESULTS 


A.  Stars 

A total  of  thirty  of  the  brighter  FK4  stars  were  observed  twice.  The 
spacing  in  declination  was  about  15°  ("from  -15°  to  +60°).  The  spacing  in 
hour  angle  was  about  2*1  (from  -4*1  to  +4*1).  Thus,  the  majority  of  the 
visible  celestial  hemisphere  was  covered. 

The  current  version  of  PEG  seeks  a distribution  of  four  stars  uniformly 
spaced  in  position  angle  relative  to  the  program  object.  If  more  than  four 
stars  are  available,  then  the  closest  or  brightest  stars  are  preferentially 
chosen.  In  addition,  the  star  closest  to  the  program  object  is  used  as  a 
control  so  that,  normally,  five  stars  are  observed.  Only  four,  however, 
are  used  in  the  modeling  discussed  in  Section  TV.  The  errors  obtained  for 
the  control  star  (e.g.,  real  position  minus  model  position)  are  then  added 
to  the  program  object’s  model  position.  If  there  are  less  than  four  stars 
(including  the  control  star)  available  within  0?75  of  the  program  object's 
position,  then  the  entire  process  is  inhibited. 

For  the  observations  of  the  stars,  the  mean  total  positional  error 
(e.g.,  mean  { [ (Aacosfi)  2 + (A6)  2]  /2 } ) was  173.  This  had  a standard  devia- 
tion of  170.  A separate  series  of  measurements  were  made  to  test  for  a 
color  or  magnitude  term.  Over  a range  of  l'l1 2 in  B-V  and  4™S  in  V no 
significant  term  appeared. 
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B.  Artificial  Satellites 


The  eight  artificial  satellites  used  to  test  the  positional  reduction 
accuracy  of  PLC  are  listed  in  Table  I along  with  their  orbital  elements. 

For  each  of  these  satellites  there  are  at  least  two  independently  determined 
sets  of  orbital  elements  available.  This  fact,  their  visibility,  and  the 
large  fraction  of  the  celestial  sphere  they  cover  prompted  our  choice  of 
these  satellites  in  particular. 

The  optical  observations  were  interspaced  with  the  radar  observations 
to  insure  that  station-keeping  maneuvers  would  be  detected.  Tn  addition, 
since  the  radar  measurements  both  predated  and  postdated  the  optical  ones, 
the  continuity  of  the  residuals  could  be  assured.  Precise  estimates  for 
the  real  accuracy  of  PLC  (when  using  artificial  satellites  as  test  objects) 
are  extremely  difficult  to  obtain  because  of  different  models  for  the  earth's 
surface,  the  geopotential's  field,  and  small  systematic  errors  in  the  measure- 
ment of  position  by  the  radar.  Also,  our  neglect  of  planetary  aberration 
is  a -2”  systematic  effect  in  right  ascension  for  the  small  inclination 
satellites.  The  neglect  of  parallactic  refraction  is  significant  onlv  for 
LAGEOS. 

Standard  software  (IMPEST)  was  used  to  produce  orbital  element  sets  by 
using  the  radar  data  only  and  combining  all  of  the  data.  We  were  able  to 
achieve  total  positional  errors  as  small  as  5"  although  the  largest  were  IS". 
Further  examination  of  the  residuals  showed  that  whenever  the  total  positional 
error  was  larger  than  *775  it  was  almost  all  in  the  declination.  This  was 
subsequently  (and  definitely)  traced  to  a problem  with  the  declination  axis 
torque  motor.  The  problem  no  longer  exists. 
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VT.  SINGLE  STAR  CALIBRATION 

Because  telescope  time  is  valuable,  and  for  many  purposes  an  artificial 
satellite's  position  need  only  be  known  to  10-15",  we  have  developed  a 
single  star  calibration  procedure.  This  proceeds  as  PLC  does  except  that 
the  star  in  the  SAOC  closest  to  the  input  position  is  chosen  to  be  the 
reference  star.  The  program  object's  position  is  then  differentially 
corrected  for  telescope  errors  and  reduced  (currently)  to  1978.0.  This 
latter  reduction  is  done  with  an  accuracy  commensurate  with  that  in  the 
topocentric  place  itself.  A non-technical  overview  of  both  PLC  and  SSC 
is  contained  in  reference  10. 
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